Bardeen-Stephen flux flow law disobeyed in the high-T c superconductor 

Bi 2 Sr 2 CaCu20 8+( 5 



A. Pallinger 1 , B. Sas 1 , I. Pethes 1 , K. Vad 2 , F. I. B. Williams 1 - 3 , and G. Kriza 1 - 4 
1 Research Institute for Solid State Physics and Optics, PO Box 49, H-1525 Budapest, Hungary 
2 Institute of Nuclear Research, PO Box 51, H-4001 Debrecen, Hungary 
3 CEA-Saclay, Service de Physique de VEtat Condense, 
Comissariat a VEnergie Atomique, Saclay, F-91191 Gif-sur-Yvette, France and 
^Institute of Physics, Budapest University of Technology and Economics, Budafoki ut 8, H-llll Budapest, Hungary 

(Dated: March 12, 2008) 

Pulsed high current experiments in single crystals of the high-T c superconductor 
Bi2Sr2CaCu20g+^ in c-axis directed magnetic field H reveal that the a6-face resistance in the free 
flux flow regime is a solely logarithmic function of H, devoid of any power law component. Re- 
analysis of published data confirms this result and leads to empirical analytic forms for the afc-plane 
and c-axis resistivities: p a b oc H 3 ^ 4 , which does not obey the expected Bardeen-Stephen result for 
free flux flow, and p c oc H~ 3 ^ 4 log 2 H. 

PACS numbers: 74.72.Hs, 74.25. Fy, 74.25. Qt, 74.25.Sv 



Free flux flow (FFF) resistivity describes how fast the 
vortices in a type II superconductor move in the direc- 
tion of an applied force It is a measure of how the 
momentum of the superfluid is transferred to the host 
lattice via quasiparticle excitations. The velocity-force 
relation expressed by the FFF resistivity has to be taken 
into account in interpreting any vortex transport, be it 
global or local. Since the primary source of dissipation 
in a type II superconductor in magnetic field is vortex 
motion, FFF resistivity is also of great importance for 
technical applications. 

A transport current exerts a Lorentz-Magnus force on 
a vortex and if other forces like vortex-defect interaction 
(pinning) are negligible (i.e., the vortex motion is "free" ), 
the velocity-force relation can be inferred from the resis- 
tivity pfff- The Bardeen-Stephen (BS) law Q states 
that pfff is proportional to the density of vortices and 
therefore to the magnetic field H: 



Pfff 



7Pn (H/H c2 f, (3=1, 



(1) 



where p n is the normal state resistivity, H C 2 the upper 
critical field, and 7 a constant w 1. This law has been 
experimentally established Q for a number of conven- 
tional superconductors. In high-T" c materials the quasi- 
two-dimensional (2d) electronic structure, the nodes of 
the d-wave order parameter, and the structure of the vor- 
tex system may potentially influence FFF. The motion 
of 2d "pancake" vortices along their well conducting ab 
plane leads to dissipation by flux flow, whereas in the 
poorly conducting c direction dissipation is governed not 
by flux flow but by tunneling between weakly coupled ab 
planes. A quasiclassical calculation Q suggests that the 
BS law is valid also in the d-wave case. 

Experiments to test the validity of Eq. JT]) in high-T c 
superconductors — and especially in Bi2Sr2CaCu20s+5 
(BSCCO), the model system of this study — are contra- 
dictory. Data on low frequency transport in single crys- 



tals [f| and thin films 0, [7 1 as well as microwave and 
millimeter wave impedance 8| are inconsistent with one 
another, agreeing only that the BS law is not obeyed. The 
resistivity often resembles a sublinear power law in field. 
The situation is similar in other high-T c materials with 
the notable exception of the results of Kunchur et al. Q 
who find agreement with BS law in thin film resistivity 
measurements in YBa2Cu3C>7. However, the resistance 
they measure does not saturate, i.e., it increases with 
current, up to the highest current they use, leading to 
an uncertainty in the value of pfff ■ A re-analysis of the 
differential resistance indicates again a sublinear field de- 
pendence of pfff- 

The main difficulty in measuring the velocity-force re- 
lation is to take account of the pinning force about which 
one has little detail. One way is to model pinning to 
interpret the surface impedance arising from local vor- 
tex motion. Another approach is to create experimental 
conditions where pinning is irrelevant as occurs in a true 
(unpinned) vortex liquid. We extend our experiments to 
the non-ohmic regime by applying sufficiently high cur- 
rent that the pinning force is negligible compared with 
the Lorentz-Magnus force from the transport current. 

To this end we have made pulsed high-current trans- 
port measurements on single crystal BSCCO with elec- 
trode contacts on the face parallel to the well conduct- 
ing ab planes in a c-axis directed magnetic field. The 
global single crystal resistance measured on the afo-face 
in the ohmic regime and the asymptotic high-current dif- 
ferential resistance in the non-ohmic regime show the 
same logarithmic magnetic field dependence devoid of 
any power law. We combine this result with published 
data from other experiments 1(J 11 1 and set up empiri- 
cal functional forms for the local resistivities p a b and p c , 
valid over a broad range of temperature and field in the 
vortex liquid phase. The most striking and important of 
our conclusions is that for BSCCO (3 — 3/4 in Eq. (1). 
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FIG. 1: Typical voltage-current characteristics at selected 
temperatures in B — 3 T. Inset: current dependence of dif- 
ferential resistance dV/dl at the same temperatures. At high 
currents dV/dl saturates at R a b- 



We selected for experiment three single crystals from 
3 different batches of Bi2Sr2CaCu20s+5 with typical di- 
mensions 1 x 0.5 x 0.003 mm 3 , the shortest corresponding 
to the poorly conducting c axis. All were close to optimal 
doping with a resistance-determined critical temperature 
T c w 89 K and transition width about 2 K in zero field; 
the diamagnetism in a 1 mT field set in progressively 
below T c to near 100% at low temperature. 

Voltage-current (V-/) response was measured in the 
usual four-point configuration on an ab face in perpendic- 
ular magnetic field with two current contacts across the 
width near the ends and two point voltage contacts near 
each edge of the same face. The contacts were made by 
bonding 25 fxxsx gold wires with silver epoxy fired at 900 K 
in an oxygen atmosphere resulting in current contact re- 
sistances of less than 3 Q. To avoid significant Joule 
heating, we employed short (< 50 /is) current pulses of 
isosceles triangular shape at 0.2 to 1 s intervals. Techni- 
cal details and the issue of Joule heating are treated in 
Ref. with the conclusion that the temperature change 
in the area between the voltage contacts is negligible for 
the duration of the pulse. 

Typical V-I characteristics at different temperatures 
in a field of 3 T are shown in Fig. [T] Above a temper- 
ature Ti in < T c the V-I curves are linear (see the lower 
inset of Fig. for the field dependence of the characteris- 
tic temperatures). Below Ti; n nonlinearity develops and 
the I — > resistance decreases faster than exponentially 
until it becomes unmeasurably small even with the most 
sensitive technique. At low temperature dissipation sets 
in abruptly at a threshold current 7th! for higher tem- 
peratures a marked upturn in the V-I curve (a "knee") 
is seen at 7k i$ Ith- Throughout the nonlinear range the 
differential resistance increases with increasing current; 
for currents several times ith or Jk it saturates (becomes 
current independent) at a value dV/dl — R a b as shown 
in the inset of Fig. [TJ Since ith and Ik are hallmarks 
of depinning, the current-independent R a t observed at 
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FIG. 2: High-current differential resistance R a t, at several 
temperatures as a function of the logarithm of magnetic field 
normalized to the upper critical field H C 2(T). Upper inset: 
Same data on a linear field scale. Lower inset: Phase diagram 
measured on same samples except for the first order transition 
Tfot taken from Ref. [14J for comparison. Ti rr refers to mag- 
netic irreversibility, T2nd to second magnetization peak, Tun 
to the beginning of linear V-I. H c2 = 120[1 - (T/T c ) 2 ] tesla. 
The vertical black (grey) lines show the range of full (open) 
symbols in the main panel. The open circle is crossover in the 
field dependence of R a b- 



many times this current suggests that at these high cur- 
rents pinning is irrelevant and R a b reflects FFF [13j . The 
focus of this article is the behavior of R a b- 

The upper inset of Fig. [5] shows the field dependence of 
R a b for the temperatures indicated in the phase diagram 
of the lower inset. The resistance is field and tempera- 
ture independent at low temperature. With increasing 
temperature there is a crossover to a field-dependent be- 
havior at T co (open circles in the phase diagram) situated 
between the magnetic irreversibility line and Tn n . Hav- 
ing in mind that in the BS law the characteristic field 
is H c2 , we interpolate the upper critical field using the 
form H c2 = (120 tesla)]l - (T/T c ) 2 ] constructed from 
dH c2 /dT Tc = -2.7 T/K [15] and use it to plot R ab against 
H/H C 2 on a logarithmic field scale in the main panel of 
Fig. O The high-temperature curves all collapse into a 
master curve representing a logarithmic field dependence: 

R ab (H, T) = R n ab [l + a \og{H/H c2 {T))] , (2) 

where i?" b is the zero field normal resistance at T c and 
a is a constant. This scaling only contains the tempera- 
ture through H C 2(T). Equation ^ provides an excellent 
description of all three samples; for the parameter a we 
find 0.16, 0.19 and 0.21, essentially the same values, in- 
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sensitive to the presumably different disorder between 
samples. 

We emphasize that R a b cannot be compared directly 
to the BS law because the strong anisotropy of the elec- 
tronic properties makes the current distribution very in- 
homogeneous and R a b reflects both a6-plane and c-axis 
properties. However, if the sample is thick in the c direc- 
tion, a simple scaling argument for a current independent 
local resistivity tensor yields R a b = A^J p a bp c where A is 
a geometrical factor. This relation is valid in the linear 
region T > Tu n (B) but also below Tn n if the current is 
sufficiently high that the current density is well in the up- 
per differentially linear portion of the response near the 
top surface of the crystal [3] • The analysis and experi- 
mental checks of Ref. [l(| indicate that with the sample 
size, shape, and contact geometry used in their and our 
single crystal ab plane studies in BSCCO, the thick sam- 
ple limit provides a good description. 

Independent confirmation of our results emerges from 
analysis of other experiments. Although high-current 
data are absent, the fact that the V-I curves are linear 
for T > Tu n (H) allows comparison with low-current data 
in this temperature and field range. In Fig.[3ja) we show 
R a b vs. H curves extracted from the R a b vs. T data taken 
by Busch et al. [l(| in different magnetic fields. Excellent 
agreement with Eq. {2J is seen for T > Ty ln (full symbols 
in the figure) with a = 0.23. 

Busch et al. (Io| were able to disentangle p a b and p c by 
using data from two additional contacts on the bottom of 
the crystal. p a b and p c results extracted from their work 
for a series of temperatures are shown as a function of 
H/H C 2 in Fig.^b) and (c). In the temperature and field 
range T > Tu n (H) both quantities individually exhibit 
H/H C 2 scaling. The in-plane resistivity does not agree 
with the (3=1 BS law of Eq. (1), but is well described 
by a/3 = 3/4 exponent (best fit (3 = 0.75 ±0.02). Having 
definite analytic forms for both R a b and p a b, we can use 
the relation (R a b/A) 2 = p a bPc to write an expression for 
the c-axis resistivity. In summary: 



Pab = Pfff = p2 b (H/H c2 f, 
p c = ^(i//i? c2 )- /3 [l + alog(iJ/i? c2 )] 2 , 
P = 3/4, a = 0.2 . 



(3) 
(4) 



The prefactors p™ b and p™ are in good agreement with 
the respective normal resistivities at T c . 

Are these forms corroborated by other types of mea- 
surement? In principle p a b can be measured in thin films 
where the current density is expected to be homoge- 
neous. In Fig. G3h) we show the 77 K thin film resis- 
tivity obtained by digitizing the V-I curves in Ref. 
Data above about 1 T are reasonably well described by 
a H 3 / 4 dependence, but a closer look reveals that the 
log Pab vs. log H curves are concave from below at every 
H, i.e., there is a systematic deviation from power law. 
This "logarithm like" (but not logarithmic) dependence 
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FIG. 3: Single crystal resistance R a t, (a); ab plane resistivity 
Pab (b); and c-axis resistivity p c (c) as a ruction of magnetic 
field normalized to the upper critical field H C 2- Data from 
Ref. For full symbols T > Ti in . Solid lines in panels (a) 
to (c) are fits to Eqs. (2) to (4), respectively. In panel (b) 
77-K thin film data from Ref. [TJ is also shown by triangles. 
Inset of panel (c): c-axis conductivity data from Ref. [ill ]. 
tj c = 1/pc and ao a constant. The solid line is a fit to Eq. (3). 



is shared with other thin film results [6( but there are sig- 
nificant quantitative differences between data measured 
by different groups. A possible reason is that macro- 
scopic defects like steps on the surface or mosaic bound- 
aries force c-axis currents and the measured resistance is 
a sample-dependent combination of p a b and p c . 

The expression for p c reproduces well the maximum (at 
ff max = H c2 exp(8/3 — 1/ce) ~ 0.1i? C 2) observed in the 
high field c-axis magnetoresistance 
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and the overall 
field dependence of these independently measured data 
is very well described by Eq. (|4]). We demostrate this in 
Fig. [3](c) where we plot 70 K data for a c (B) — ao from 
Fig. 5 of Ref. ll| where a c — l/p c and the constant ao 
is interpreted as the zero-field quasiparticle conductivity. 
Using our estimate of H c2 (T) » 46 T for T = 70 K, 
Eq. ([4]) fits the measured data with parameters <r" = 
1/p™ = 6.6 (kfWT 1 , ao = 3.6 (MWT 1 and a = 0.20 
to obtain the curve indicated by the continuous line in the 
figure. The value of a is in excellent agreement with that 
inferred from R a b measurements. It should be pointed 
out, however, that ctq is significantly smaller than the 
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value do ~ 8 (kficm) -1 inferred in Rcf. llj. In terms of 
resistivities, this means that p c decreases more slowly in 
high fields than described by Eq. ^ with (3 = 3/4 and 
is in fact best described with an exponent (3 = 0.51. 

In the high-current limit the same form for R a b(H) also 
holds below T\[ n (H) where the V-I curves are nonlinear. 
Since no change in the behavior of R a b is observed when 
the T\i u {H) line is crossed, it is reasonable to assume 
the same for p a b and p c . In the low field direction a 
lower limit of the validity of Eq. ([2]) is the zero of the 
equation at H /H c2 = e~ x l a ~ 10~ 3 — 10~ 2 , higher but in 
the order of the first order transition in the static vortex 
system. In the high-field direction Eq. ([2]) is valid up to 
the highest field « 0.3H C 2 we investigated. 

The temperature T co of the crossover from R a b — const 
to R a b oc log H is distinctly higher than the onset of mag- 
netic irreversibility at T; rr , and also above the vortex glass 
transition T g « T; rr inferred from scaling analysis [la ] of 
the V-I curves. On the other hand, no change in the 
behavior of R a b is observed when the T- m (H) and T g (H) 
lines are crossed. This suggests that because the pin- 
ning potential is smoothed at high velocities, the phase 
diagram |17| of the far-from-equilibrium dynamic vor- 
tex system [18[ is different from that of the unperturbed 
thermodynamic phases. Since R a b behaves the same in 
the pinned (T < Ti; n ) and unpinned (T > Ty ul ) liquid 
phases, we propose that the unpinned phase, otherwise 
observed only above Ty ln , may be restored in the range 
T co < T < T\i n . Then T co may approximate the melting 
transition in a hypothetical defect-free crystal. 

Our most robust finding, invariably observed not only 
in our 3 batches but also in the data of Ref . dl| , is the 
logarithmic field dependence of the high-current single 
crystal resistance R a b- Although a power of H factor is 
expected both in p a b [l| and p c [19(, no such factor is 
present in R a b oc p a bp c . The most likely reason is that 
the power-law factors in p a b and p c cancel (exponents 
3/4 and -3/4 in our analysis). The cancellation is very 
accurate; we estimate that a power law factor with ex- 
ponent as low as 0.1 could be observed in our R a b data. 
Moreover, because we find no logarithmic correction to 
Pab, the logarithmic dependence of R a b is carried by p c . 

Arguing that both p a b and a c are proportional to the 
quasiparticle density of states at the Fermi level, N(0), 
it cancels in the product p a bPc- In conventional super- 
conductors N(0) is proportional to the number of vor- 
tices therefore to H, leading to the ii-linear resistiv- 
ity of the BS law. In nodal gap superconductors near- 
nodal quasiparticles lead to a sublinear dependence; for 
line nodes N(0) oc H 1 / 2 [2(J, as evidenced in recent 
low-tempcrature thermodynamic measurements 21| . Al- 



ent structure factors of the solid and liquid phases. 

Nonlocal effects 22, 23| may influence the evaluation 
of the 6-contact measurements [l(| and therefore the va- 
lidity of Eqs. ([3]) and f4} (but not of Eq. ([21)). This seems, 
however, unlikely in the light of the goo d agreement of 
p c inferred from independent a6-plane |l0| and c-axis [ll| 
measurements and of the broad temperature range of va- 
lidity of Eq. ©. 

In conclusion, we have set up empirical rules for the 
analytic form of single crystal resistance as well as for 
the ab plane and c axis resistivities in the high-current 
free flux flow limit in the vortex liquid state of BSCCO, 
valid over a broad range of temperature and field. Both 
the logarithmic field dependence of the single crystal re- 
sistance and the 3/4-power law in the a6-plane free flux 
flow resistance are in disagreement with the current the- 
oretical understanding of high-T c superconductors. 
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EPFL laboratory in Lausanne have contributed in a very 
essential way to sample preparation and characterization. 
Finally we acknowledge with gratitude the Hungarian 
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though delocalized near-nodal quasiparticles are not ex- 
pected to contribute significantly to pfff because of the 
weak spectral flow force [3] they experience, the result 
may be different in the diffusive limit in the liquid phase. 
A possible reason for the (3 = 3/4 exponent is the differ- 
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